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Some recent work on the thermodynamic behavior of the matrix 
model of M-theory on a pp-wave background is reviewed. We examine 
a weak coupling limit where computations can be done explicitly. In 
the large N limit, we find a phase transition between two distinct 
phases which resembles a "confinement-deconfinement" transition in 
gauge theory and which we speculate must be related to a geometric 
transition in M-theory. We review arguments that the phase transition 
is also related to the Hagedorn transition of little string theory in a 
certain limit of the 5-brane geometry. 

1 Prologue 

Ian Kogan was a great friend and I will miss him dearly. Part of his 
journey from the Soviet Union to Oxford passed through Vancouver 
where he spent a few years. I have great memories of that time. 

Among Ian's very broad range of scientific interests was a contin- 
uing fascination with critical behavior of string theory at high tem- 
perature. One characteristic of strings is that their density of states 
increases exponentially at large energies, 



The constant Tjj is called the Hagedorn temperature. A consequence 
of this large density of states is that, depending on the exponent, a, 
string theory has either an upper limiting temperature or a phase 
transition. 
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Perhaps my favorite of all of Ian's scientific works is an old result [1], 
(found independently in [2]) about an interpretation of the Hagedorn 
temperature in string theory. In that work he noted an analogy be- 
tween the Hagedorn behavior of strings and the Kosterlitz-Thouless 
phase transition for the unbinding of vortices in the world-sheet sigma 
model. He interpreted the latter in string theory as a disintegration 
of the worldsheet by condensation of vortices. Characteristic of Ian's 
work, this very original and fascinating idea seemed well ahead of its 
time, In all likelihood, its full import has yet to be realized. 

Ian's interest in the high temperature behavior of strings continued 
throughout his career. Some of his recent work explored the use of 
the AdS/CFT correspondence to understand the phase structure of 
string theory at high temperatures [3, 4, 5]. 

In this Paper, which I dedicate to Ian, I will discuss some of my 
own recent work on similar topics. 



2 Motivation 

The basic degrees of freedom of string theory and M-theory are thought 
to be known and encoded in the BFSS [6] matrix model. The model 
is super symmetric matrix quantum mechanics with action 



/ 



dtTv 



(2) 

where i, j = 1, ...,9 and all degrees of freedom are N x N Hermitian 
matrices. This is a gauge invariant theory with covariant time deriva- 
tive D = 4: — i[A, .... The gauge theory coupling constant is given in 
terms of the null compactification radius R and the eleven dimensional 
Planck length i v \ (or the ten dimensional IIA string coupling g s and 
string length l s = \fa!) by 



9s£ 3 s , i P i = gl /3 ls , R = g s £s (3) 



This model has three uses. It is conjectured to describe a discrete 
light-cone quantization of M-theory [6] where R is the compactification 
radius of the light-cone, there are N units of light-cone momentum and 
l p is the Planck length of 11-dimensional supergravity. Secondly, and 
historically a little earlier [7], with parameters suitably re-identified, it 
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describes the low energy dynamics of a collection of N DO-branes of 
type IIA superstring theory. Finally, it is a matrix regularization of 
the light-cone action for the 11-dimensional supermembrane[8]. 

One motivation for understanding the behavior of matrix models 
such as the BFSS model at finite temperature comes from the con- 
jecture that their finite temperature states are related to black hole 
states of type IIA super gravity [9, 10]. This idea was studied in a series 
of papers by Kabat and Lowe[ll]. They begin with the Beckenstein- 
Hawking entropy of a black DO-brane solution of IIA supergravity - 
the area of its even horizon in Planck units. They convert the entropy 
to the free energy, which they then write in terms of gauge theory 
parameters to obtain 



A derivation of (4) from the matrix model would be an important 
result, a first principles computation of non-extremal black hole en- 
tropy using string theory. However, one would expect to find (4) in 
a low temperature and therefore strong coupling limit of the matrix 
model, making it inaccessible to perturbation theory. A variational 
technique was applied and claimed approximate agreement with the 
formula over some range of temperature [11]. 

The formula (4) is remarkable in three respects. First, it has the 
correct dependence on N and the 'tHooft coupling gy^N to be the 
leading order of the 'tHooft limit of the gauge theory. If it could 
be derived in a perturbative expansion, it would obtain contributions 
only from planar Feynman diagrams. This means that the 'tHooft 
limit should be part of the limiting process that would extract the 
classical physics of black holes from the matrix model. 

Secondly, the scaling with N in (4) is as if this 0-dimensional gauge 
theory were in a de-confined phase. This is particularly true if we as- 
sume that we are taking the 'tHooft limit. The use of the word "decon- 
fined" in a theory where there is no spatial extent over which particles 
can be separated must be justified carefully. The gauge theory has a 
Gauss law constraint so that quantum states of the Hamiltonian must 
be singlets under the gauge symmetry. The gauge field in (2) enforces 
this constraint. The number of singlets at a given energy do not scale 
like N 2 , rather they are of order one. As an example of this, consider 
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a matrix harmonic oscillator with Hamiltonian 

d 

tf = 5>Tr(aia;) (5) 

i=i 

and matrix-valued creation and annihilation operators with algebra 

= S^S ad 6 bc (6) 

States are created by ot£ b operating on a vacuum |0 >. To get a state 
with energy E = nu; we must act with n creation operators. 

The analog of gauge invariance is to require a physical state con- 
dition of invariance under the unitary transformation 

— ► uaftvt 

where 

u e U{N)/U{1) 

We assume that the vacuum state is invariant under this gauge trans- 
form. Then, physical states are created by operating with invariant 
combinations of creation operators. In the limit N — > oo all such 
combinations are traces 



where the energy is 

E = oj(m + 2n 2 + 3n 3 + . . .) 

The number of these traces with a fixed energy, E, does not scale like 
iV 2 as N — > oo, instead it approaches a constant as N is taken large. 
Thus for normal thermodynamic states, one would not expect the free 
energy to be of order N 2 . 

However, the number of independent traces does increase rapidly 
with E. It has been shown[12] that, in the large N limit, the oscillator 
has a Hagedorn-like density of states at high energy, 

p(E) ~ 

where the Hagedorn temperature is 

T H = uj/lnd 



Tr (a 



31 32 
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"3 



|0 > (7) 
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A similar result has been found for weakly coupled Yang- Mills theory[13, 



At temperatures higher than the thermodynamic canonical 
ensemble does not exist. It could be made to exist by keeping N 
large but finite. That would cut off the exponential growth in the 
asymptotic density of states at some large energy. Then we could 
consider a temperature that is greater than Tjj- Both the energy and 
entropy would be dominated by states at and above the cutoff scale. 
Then, the divergence of the free energy ~ N 2 occurs as we take the 
limit N — > oo at constant temperature (noting that the Hagedorn 
temperature does not depend on N). 

As we shall show in the following, a behavior like this can indeed 
by found in the matrix model. In more conventional terms, it occurs 
as a large N Gross- Witten type of phase transition [16, 17] which is 
familiar in unitary matrix models. It is this behavior that we call 
"deconfinement" . At this deconfinement transition, limjv-xxi F/N 2 
jumps from being zero to of order one. 

In an adjoint gauge theory such as (2), there is an order parame- 
ter for confinement, the Polyakov loop[18, 19]. It is the trace of the 
holonomy of the gauge field around the finite temperature Euclidean 
time circle, 



This operator gets a nonzero expectation value when a gauge theory 
is deconfined. An interesting question is whether it has a nonzero 
expectation value in the BFSS matrix model. In such a low dimen- 
sional theory, it can only have a nonzero expectation value when iV is 
infinite. Indications from weak coupling computations [20] are that it 
has. 

The third remarkable fact about (4) is, though this formula is 
thought to apply to the black hole only for a range of temperatures [21], 
the expression at low temperature is reminiscent of critical scaling with 
a critical temperature T = and a simple, rational critical exponent. 

In spite of the simplicity of these interesting features, there is no 
analytic derivation of the formula (4) from the matrix model. One 
of the difficulties in finding a derivation is the intractability of the 
model itself. These difficulties are well-known from previous attempts 
to analyze its thermodynamics [20]. 

Before we continue with the matrix model, we comment that, if we 
analyze the thermodynamics of M-theory in the rest frame, we would 



14, 15]. 




(8) 
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form the partition function with a Boltzman distribution using the 
energy E = j-(P+ + P~), 

£ = £ e -iV/^ e -F/V2T (g) 
N 

where the matrix model free energy is defined by 

e -F/T = ^-H/T (1Q) 

In (9) we have traced over the eigenstates of the light-cone momentum. 
This gives the sum over matrix model partition functions for each N 
with the exponential factor e~ N ' RT . The convergence of the sum is 
clearly dependent on the nature of the large TV" limit. If there is a 
sector of the matrix model where this limit is like (4), because of the 
negative sign, the sum over ./V diverges. It is tempting to associate 
this with the non-existence of thermodynamics of a theory of quantum 
gravity on asymptotically flat space - because of the Jeans instability 
the space is unstable to collapse to black holes. If there were a phase 
where the free energy did not become negative and with magnitude 
growing faster than N, it would be a stable phase. 

There are several known behaviors of matrix models in the large N 
limit. For example, there is the Dijkgraaf-Verlinde-Verlinde[22] limit 
of matrix string theory. That is a strong coupling limit which kills 
the off-diagonal degrees of freedom of the matrices. The remaining, 
diagonal degrees of freedom are N in number and it can be shown 
explicitly [23, 24, 25] that the free energy is negative and is proportional 
to N. Of course, this is just the correct behavior for a string theory, 
there will be a Hagedorn temperature where the large N terms in 
the sum in (9) go from being exponentially suppressed to growing 
exponentially. There are other versions of matrix string theory[26, 27, 
28, 29] based on two dimensional Yang-Mills theory where one would 
expect a similar behavior. 

Another limit where we have a quantitative estimate of the large 
N behavior is the 'tHooft limit where N is taken to infinity at the 
same time as qym is taken to zero. Technically this would be done by 
re-defining A = <7y M iV and holding A fixed as when we sum over N 
in (9). Then, the phase transition that we discuss here is somewhat 
more violent than the Hagedorn behavior in string theory. The matrix 
part of the free energy at large N changes from a negative constant 
to a negative constant times N 2 . The linear in N exponent of the 
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momentum part doesn't compete with the iV 2 -growth of the matrix 
model contribution. 



3 PP-Wave matrix model 

Recently a variant of the matrix model which is conjectured to de- 
scribe a discrete light cone quantization of M-theory on a pp-wave 
background has been formulated [30]. It is a 1-parameter deformation 
of (2), 



-/ 



S = / dtTr 



~mi {xa)2 ~ 4r {x ' 1 ' )2 ~ 4^ " ^ ab ° xaxbxc 



(11) 



where the indices a, b, .. = 1, 2, 3 and i! = 4, 9. This matrix model 
reduces to the BFSS model if we put the parameter \i to zero and can 
be considered a one-parameter deformation of it. The main difference 
between the two is that the action in (2) has (super) symmetries iden- 
tical to those of the residual invariance of 11-dimensional Minkowski 
space in light-cone quantization whereas (11) has symmetries appro- 
priate to a pp-wave spacetime. 

The matrix model in (11) has the great advantage that, unlike (2), 
it can be analyzed in perturbation theory [31]. In (2), the classical 
potential ~ — Tr ([X 1 , X^] 2 ) has flat directions, any set of matrices 
which are mutually commuting have zero energy. The behavior of the 
degrees of freedom in flat directions must be understood at the outset 
of an honest quantum mechanical treatment of the theory. In (11), 
these flat directions are removed by the mass terms. Perturbation 
theory is accurate in the limit where the mass gap \x is large. 

The pp-wave space which is a maximally supersymmetric solution 
of 11-dimensional supergravity is 

ds 2 = dx'dx* - 2dx+dx- - \^-{x a ) 2 + ^(x*') 2 J dx + dx + 

with an additional constant background 4-form flux 

-F+123 = M 
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This space is known to support a spherical membrane solution, 



x = p t ,x 



= const. , \J (x a ) 2 = 



4 



p' H 



and a spherical transverse 5-brane 



x + =p+ T , x' = const. , 0/) 2 = £ 3 ■ n ' P 



V 1 
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These objects are conjectured to make their appearance in the 
solutions of the matrix model. The membrane is found immediately 
in semiclassical quantization. The classical potential is (hereafter we 
set £ p i = 1 and measure all dimensional quantities in Planck units) 



V=*Tr 
2 



{ik** + ieabCXbXC ) 2 + \ ' Xj '0 2 + 

+ {i[x i ',x-}Y + ^\x i ') 



It is minimized by 



cl > cl 



(12) 



(13) 



Where J a is an N-dimensional representation of the SU(2) algebra, 
[J a ,J b ] = ie abc J c . In addition, the classical solution for gauge field 
must obey the equation 

If J a is an irreducible representation of SU(2), by Schur's Lemma, 
A c \ = 0. The gauge symmetry is realized by the Higgs mechanism. 
When the representation is reducible, there are gauge fields which 
commute with the condensate. This part of A c \ remains undetermined 
and must still be integrated over, even to obtain the leading order in 
the semi-classical approximation to the partition function. 

The configurations in (13) are fuzzy spheres, which are matrix 
regularizations of the membranes. The 5-branes on the other hand do 
not seem to appear in the perturbative states of the matrix model. It 
has been conjectured [32] that the 5-branes indeed appear as the large 
N limit is taken in a certain way. 
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First of all, J a need not be an irreducible representation, but can 
contain a number of irreducible components, 



J a = 



/sf \ 

sfj 

si 

V 8 % J 



To get a (multi-)membrane state, the large N limit is taken by holding 
the number of representations fixed and sending the dimension of each 
of the representations to infinity. 

To get a five-brane state on the other hand, we hold the dimensions 
of the representations fixed and repeat them an infinite number of 
times to get the large N limit. 

An important difference between these limits is in the realization 
of the gauge symmetry. In the classical sectors, the gauge symmetry is 
partially realized by the Higgs mechanism, with the residual symmetry 
being that which interchanges representations of the same size. In a 
membrane state the residual gauge group thus has finite rank, whereas 
in a 5-brane state its rank always goes to infinity. The single 5-brane 
state is X a = whereas the state with k coincident 5-branes has the k- 
dimensional representation repeated N/k — > oo times. A state with k 
non-coincident 5-branes has largest representation k-dimensional and 
a number of smaller representations all repeated an infinite number of 
times. 

The effective coupling constant which governs a semi-classical ex- 
pansion about one of the classical ground states is 



where n is the rank of the residual gauge group. The 5-brane limit 
is where we are required to take the weakest coupling limit. It is 
the conjecture of ref.[32] that the membranes and all other degrees of 
freedom decouple in this limit and it isolates the internal dynamics of 
the 5-brane. 



Now, let us consider a perturbative expansion of the pp-wave ma- 
trix model at finite temperature. The partition function is the path 




4 Perturbative expansion 
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integral 

Z = J [dA][dX l \[d4>]e- drL[A ' Xl ^ ] 

where L is the lagrangian with Euclidean time and /3 = 1/T is the in- 
verse temperature. The bosonic and fermionic variables have periodic 
and antiperiodic boundary conditions, respectively 

A(t + /3)=A(t) , X i (r + P)=X i (r) , ^(r + /?) = -^(r) 

Since the fermions are antiperiodic, these boundary conditions break 
supersymmetry. Of course this is expected at finite temperature. 

We begin by fixing the gauge. It is most convenient to make the 
gauge field static and diagonal, 

~i~A a i, = , A a i) = A a 5 a i) 
dr 

The remaining degrees of freedom of the gauge field are just the time- 
independent eigenvalues A a . 

The Faddeev-Popov determinant for the first of these gauge fixings 
is this gauge fixing is 

det '(-|(-^ +i ^-^)) < 14 > 

where the boundary conditions are periodic with period (3. The prime 
means that the zero mode of time derivative operating on periodic 
functions is omitted from the determinant. The Faddeev-Popov de- 
terminant for diagonalizing the gauge field is the familiar vandermonde 
determinant, 

n \ A a - M 

Together the two determinants are 




where the prime on the determinant indicates that the static mode is 
omitted. 1 

Using zeta-function regularization, 

det '(-!H 
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If we expand about the classical vacuum corresponding to the sin- 
gle 5-brane, X"j = = X l cl , we find the partition function in the 1-loop 
approximation is 



J dA a \{ 



det' (-d/dr) det (-D ab ) det 8 (-D ab + \ ) 
J£ det 3 / 2 (-DI + £) det 3 (-^ b + g) 



(15) 



where -D a ^ = ^ — i(Ai — Using the formula 

/ d \ , /?w 

det — - — hw = 2 sinn — 
V dr J 2 

with periodic boundary conditions and 

, , d \ Blo 
det — - — h to = 2 cosh ■ 



dr J 2 
with antiperiodic boundary conditions, we can write 



1 r 1 / 



1/2 V 2vr / 11 [1 — e - J 9M/3+*/9(A a - J 4 6 )]3[i _ e - / 9/x/6+i/3(A a -A 6 )]6 

(16) 

The factor of 1/N\ is the volume of the residual discrete gauge group 
which permutes the eigenvalues. When N is finite, it might be possible 
to do this integral using the method of residues. 

However, to apply to the 5-brane, we require the integral when 
N — > 00. There are N integration variables A a and the action, which 
is the logarithm of the integrand is generically of order TV 2 which is 
large in the large N limit. For this reason, the integral can be done by 
saddle point integration. This amounts to finding the configuration of 
the variables A a which minimize the effective action: 

S cS = E (" Ml " e lf3 ^- A ^} - 81n[l + e -/W4+tf(A,-A6)] + 

a^b 

+ 3 ln[l - e ~M3+iP(Aa-A b )^ + 6 ln [! _ e -/3»/6+if3(A a -A b )^ (17) 

To study the behavior, it is illuminating to Taylor expand the loga- 
rithms in the phases (this requires some assumptions of convergence 
for the first log) 

Ses = y »l ) e 6e be | Tre ^A|2 (lg) 

^1 
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Each term contains the modulus squared of a multiply would Polyakov 
loop (8). When a coefficient becomes negative, the loop condenses. In 
fact, as we raise the temperature from zero (and lower (3 from infinity), 
the first mode to condense is n = 1. This occurs when 

T c = ^ ln3 ~ -0758533^ 

The condensate breaks a symmetry under changing the phase of the 
loop operator. 



5 A closer look at the phase transition 

At temperatures greater than Tc the eigenvalues A a distribute them- 
selves so that they are clustered near a particular point on the unit 
circle. To examine the possibility, we consider the equation of motion 
for the eigenvalues, 

uj(e 0+ z) + oj(e°-z) + 8u(-r 3 z) + 8w(-r~ 3 z) = 

= 3u(r 4 z) + 3to(r- A z) + 6uj(r 2 z) + 6w(r~ 2 z) (19) 

where r = eP^I 12 and the resolvent is defined as 

u( z ) = ± V -o A (20) 

a=l 

uj{z) is holomorphic for z away from the unit circle and has asymptotic 
behavior, uj(oo) = 1 and co(0) = —1. In the large N limit the poles 
in u>(z), which occur at the location of the elements a , coalesce to 
form a cut singularity on a part or perhaps all of the unit circle. lo(z) 
remains holomorphic elsewhere in the complex plane. 

We must remember that equation (19) is valid only when z is 
one of the gauge field elements e l ^ Aa . In that case, the sum in (20), 
which turns into an integral in the large N limit, must be defined as a 
principal value. In (19) this is gotten by averaging over approaching 
the unit circle from the inside and from the outside. 

It is easy to find one exact solution of (19). If we consider the case 
where a are uniformly distributed over the unit circle, so that the 
sum in (20) is symmetric under z — ► e ie z we can average over the 
symmetry orbit to get 

!:!<! pD 
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The result is certianly holomorphic everywhere away from the unit 
circle and is discontinuous on the entire unit circle. 

The resolvent (21) is always a solution of (19) for any value of r. 
This is the symmetric, confining solution of the matrix model, where 
the Polyakov loop operator, whose expectation value is a particular 
moment of lo(z) for large z, vanishes. We would expect that this 
confining solution is only stable if the temperature is low enough. At 
some critical temperature it becomes an unstable solution and there 
should be other solutions which have lower free energy. 

The confining phase which we discuss above is stable when r > 3 
or r < 1 /3. When r = 3 or r = 1 /3, we can find a 1-parameter family 
of solutions, 

/ \ ( -1-az \z\ < 1 . . 

= i ! , / i ' i ' -, (22) 



1 + a/z \z\ > 1 

This is an acceptable solution when \a\ < 2. 2 If we plug it into eqn.(19) 
and assume that r > 1, we obtain 

(Sr- 4 + 8r~ 3 + 6r~ 2 -l)(z-l/z) =0 (23) 

which is solved when r = 3. If we assume r < 1 we find an equation 
which is solved by r = 1/3. 

To examine the phase transition further, we expand about r = oo. 
We expect the transition to occur at r = 3 which is not really large, 
but we will see that corrections are of order 1/r 4 , at the 1-percent 
level. 

The asymptotic expansions of the resolvent is 

oo 

u(z) = l + 2'£—T h (24) 

-i z 
n=l 

oo 

^(z) = -l-2^z"r/_n (25) 

71=1 



infSAa 



where 

1 N 

a=i 

are the expectation values of the Polyakov loop operator for n wind- 
ings. If we assume that r > 1, an asymptotic expansion of the equation 



2 Hcrc, a/2 is the expectation value of the Polyakov loop operator which must be less 
than one. 
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(19) is 



^z)+u{e^z) = 2 £ (A + JL + M (%z -n _ ^n) (26) 

n=l ^ ' 

Remember that this equation is valid only when z is inside the cut 
discontinuity of u>(z) which is assumed to occur on a segment of the 
unit circle. 

In the large r limit, the right-hand-side of this equation can be 
approximated by the leading terms. It is then similar to the equations 
for the eigenvalue distributions in adjoint unitary matrix models which 
have been solved in the literature [17, 34, 35]. 

It is easy to find a solution of (26) if we truncate the right-hand-side 
by retaining only the n = 1 term. Consider the semi-circle distribution 
of Gross and Witten[16] 3 

1 {x ^ 1 'i+r\ vi+2tz+*2 



= —t{- z - z )-(TTT){ 1 + - z ) Vl + 2tz + z2 ™ 

This function has a cut singularity on the unit circle between branch 
points z± = — t ± i\J\ — t 2 where we take t in the range — 1 < t < 1. 
When t — > 1 the endpoints of the cut touch each other and the cut 
covers the whole unit circle. This is the where the Gross- Witten phase 
transition occurs in their unitary matrix model. [16] In their case, it is 



3 The spectral density is denned by 

JV 



It is normalized so that 



d0p(9) = 1 



In the large N limit, it becomes a continuous function of 8 with support on some or all of 
the interval [— tt,tt}. An example is the semicircle distribution, which is given by 



rm _/^rr+t) cos fv /2 ( 1 + *)- 4sin2 l o<sinf<yiS 

Psc\") — 



To get (27) we integrate 



Jl±l < S i n | < 1 



jo 



r z + e' 

UJ sc {z) = / d9 _ ie pac{&) 
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a third order phase transition. In the present case it is a first order 
phase transition. The solution that we found above, when r = 3, is 
just a special case of the semicircle (27) when t = 1. 

Let us explore oj sc (z) a little more. First we note that it obeys 

u sc (l/z) = -uj S c(z) 

(with the appropriate change in the sign of the square root). This 
means that the % = r/^ n = r\* n for all n. We can expand for small z, 

^ )= - 1 -^_(^- (i-'>y'+^ + ,.. m 

from which we identify 

3-t (l-tf (1 - t) 2 (5t + 1) 
?yo = 1 , Vi = ' ^ = 4 ' % = 16 ••• ( 29 ) 

We see that, there is a critical point at t = 1. At that point, ryo = 1, 
77-j-i = 1/2 and ?7|fc|>i = 0. This is precisely the value of t for which 
the edges of the cut meet, so that the cut covers the entire unit circle. 
This is also precisely the exact solution (22) which we found when 
r = 3, here with the special value a = 1/2. 

In fact, the semicircle distribution gives a good approximation to 
the solution when r is slightly less than 3. For z in the cut, 



uj sc (e 0+ z) + io sc (e°~ 



l + t\z 



1 



z 



If, for the moment, we truncate the right-hand-side of (26) to the 
term with n = 1, we see that the equation is solved by the semi-circle 
distribution when 

1 = (8r~ 3 + 3r~ 4 + 6r~ 2 ) (30) 

(l+t)(3-t) V ) y 1 

Also, remembering that t falls in the range— 1 < t < 1, we get the 
critical value of r, r cr jt = 3 by setting t = 1. (30) has a solution only 
when r < r cr jt = 3 (and, here we have assumed r > 1). 

When r = 3, the n = 2 term on the right-hand-side of (26), which 
we have ignored, contains (8r~ 3 + 3r -4 + 6r~ 2 ) = .086. Thus, we see 
that, to an accuracy of about ten percent, the semicircle distribution 
is an approximate solution of the model for temperatures just above 
the critical temperature. 
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5.1 Systematic improvement of the semicircle 

It is clear what has to be done to improve this approximation. We 
can begin with an Ansatz for the resolvent which has a single cut 
singularity placed on the unit circle 

K 

u{z) = ]T (a n (z~ n - z n ) - b n (z n - z'"' 1 ) y/l + 2tz + z 2 ) (31) 

n=l 

To get the general solution, we should consider all orders by putting 
K — > oo. An approximate solution is found by truncating at some 
order K. We will see below that this approximate solution is good 
near the phase transition. The coefficients in (31) must be arranged 
so that, in an asymptotic expansion in small z, 

1. all of the poles of order l/z K , 1/z cancel and oo(0) = —1 so, 
the asymptotic series then has the form 

co = — 1 — 2i]iz — 2r/2Z — . . . 

This gives K+l conditions that the 2K+1 parameters (a n , b n , t) 
must obey. 

2. From the above expansion, we determine the moments in terms 
of the parameters 

t) , m(an,b n ,t) , ... , r] K (a n ,b n ,t) 

3. Then we use the equation (26), with the right-hand-side trun- 
cated to order K to get K conditions 

«i = f(r)vi(a n ,b n ,t) , (32) 
02 = f{r 2 )m(.a n ,b n ,t) , (33) 
a K = f(r K )r] K (a n ,b n ,t) (34) 

where f(r) = + ^ + ^j- This gives K further equations 
which completely determine the 2K + 1 parameters of the solu- 
tion. 

For example, if we choose K = 2, requiring that the poles cancel and 
uj(0) = —1 yields the three conditions 

a 2 -b 2 = 0, a 1 -b 1 -tb 2 = 0, h(l + t) + 6 2 -(l - t 2 ) = 1 
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We can use these equations to eliminate 61,01,02 

t ~ 1 , 1 3t " 1 . 1 

bl = — b2 + —t ' ai = — h2 + —t ' 02 = 62 

Then, we can calculate the first and second moments by considering 
an asymptotic expansion of uj(z). We get 

(t + 1) 3 6 2 - 2(t - 3) 
771 = 8 

m = -Y^t-h)(i + tf + \(i-tf 

and finally, using (32) and (33), we get the equations 

^ + rii = H (t+1) V (t - 3> ) 

6 2 = /(r 2 ) (-^(3t-5)(l + t) 3 + ^(l-t)- 

Of course, we already know that these equations are solved at the 
critical point by 62 = 0, t = 1, r = 3 — > /(r) = 1. If we consider a 
value of r somewhat less than the critical value, 



In this case, 

We get 
and 



/(r) = 1 + e 



V 7 2187 162 



t = 1 - 2y/e 
187 



-e 



2000 

This demonstrates that, close to the phase transition, the semicir- 
cle distribution gives an accurate description of the de-confined phase 
and this description can be systematically corrected. It would be in- 
teresting to explore this further to determine precise thermodynamic 
properties of that phase. The next term in the series on the right- 
hand-side of (26) which we have ignored, since we truncated to order 
2, is proportional to /(r 3 ) with r ~ 3, which suggests that in the vicin- 
ity of r = 3, the error is less than one percent. However, we caution 
that this is the case only for r close to 3. When r = 2, /(r 3 ) = .11. 
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5.2 High temperature limit 



Another limit we could consider is the high temperature limit where 
r — > 1. In that case, we expect that the values of e % @ Aa are con- 
centrated near a point. In fact, the case where they are at a single 
point 



is a saddle point for all values of r. However, for r ^ 1 it has infinite 
positive energy, crossing over to infinite negative energy when r = 1. 
To see that it is a solution, we note that, in this case the eigenvalue 
support is at z = rj and therefore the variable in (19) is r/. Then 
ui{e QJr rf) +uj(e°~ri) = 0. Also oo(r s r)) + ui(-^r)) = and (19) is solved. 
It is easy to see from the action that this solution is unstable for all 
values of r except r = 1. 



5.3 Free energy 

This shows the nature of the phase transition. At the critical point, it 
will turn out that the free energy is continuous, but the expectation 
of the Polyakov loop is equal to a and is ambiguous. Just below the 
transition, the theory is approximately described by the semi-circle 
distribution for which the Polyakov loop is 1/4. So we see that it 
jumps in value from to 1/4 at the phase transition. It is for this 
reason that we expect the transition to be of first order. Indeed, by 
examining the free energy, we see that it is given by 
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= 772 E (" ln ( z « - z b) " 81n(z a - r 3 z b ) + 3ln{z a - r A z b ) + 61n(z a - r 2 z b )) 

(35) 



For the symmetric solution, where po(6) = ^ 

70 = 

When r is large, we can expand to get 



7 



N 2 



E 

a^b 



In 1 1 — z a /z h \ 



6_ + 8_ + _3 \ Za 

r 2 r 3 r 4 J z b 



(36) 



If we keep only the first term in the large r expansion (note that there 
is another term which competes with 1/r 4 which we ignore for now), 
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this is approximately an adjoint unitary matrix model of the kind 
solved in the literature [34]. It is solved by the semicircle distribution 
or the symmetric distribution. The free energy is 

7 "{^(i- v /5 ^)-/W + f + ^ n [/w( 1 + v /rT ^^ *<* 



(37) 



The value of the Polyakov loop is 
1 




jfTre i $ A ) = { 1 _(, M /T— O - " (38) 



5.4 Symmetry restoration? 

Because of the low dimensionality of the system that we are discussing, 
the symmetry breaking which occurs in the de-confined phase could 
be destroyed by quantum fluctuations. In fact, it would generally be 
the case in theories with local interactions. 

For example, if N is finite, symmetry breaking is not possible. The 
phase transition that we have discussed here can only occur when N 
is infinite. Mathematically, we can think of large N as the analog of a 
large volume limit in a statistical mechanical system. If N is large but 
not infinite, the symmetry is not broken in a mathematical sense but 
the decay rate of a non-symmetric state is exponentially suppressed 
in the volume, in this case ~ e~ N "\ 

The deconfined solution has a spectral density p{6). Because of 
symmetry of the problem under replacing 9 by #+constant, there 
would be a zero mode of the linear equation for the fluctuations of 
p(9), with wave-function ip(0) ~ 4sp(0)- This mode would provide 
the motion which would restore the symmetry. However, in the case 
of the semi-circle distribution, because of the square-root singularity 
at the edge of the distribution, this function is not square-integrable, 
and therefore not normalizable. 



6 A stack of 5-branes 

The fluctuation spectrum of the matrices about a stack of k 5-branes 
is known. If, rather than the trivial vacuum, we had chosen the 
one where the fc-dimensional representation of SU(2) is repeated N/k 
times, the residual gauge invariance would be SU{N/k). The gauge 
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field would have a classical solution with k x k unit matrices times 
gauge fields A a representing the whole blocks. The spectrum is known [32, 
33] and we can again get an estimate of the Hagedorn temperature[12] 

T H (k = 00) = .073815...// 

We see that the temperature is reduced only slightly in this case. 
We take this as evidence that the Hagedorn temperature in this weak 
coupling limit is insensitive to th enumber of 5-branes. This seemingly 
contradicts the /c-dependence of the Hagedorn phase transition of little 
string theory which has been computed using holography and which 
behaves like Tc ~ \fk. 

An explanation for this contradiction can be found in the limit that 
we are using. It is a large N 'tHooft limit and a further expansion 
in weak 'tHooft coupling. In this limit, if we were to translate the 
parameters of our model to those that would describe the NS five- 
brane in type II string theory, the radius of the spherical five-brane 
would be 

r 2 

a' 

This means that we are expanding about a small, highly curved five- 
brane, whereas the usual holographic result for the Hagedorn temper- 
ature is for a large flat 5-brane. This is a similar difficulty as the one 
which appears in the AdS/CFT correspondence in general. There, 
the analog of the matrix model, which is maximally supersymmetric 
Yang-Mills theory, can be readily analyzed only in the weak coupling 
limit using an asymptotic expansion in A. Super gravity and holog- 
raphy give tools which compute the strong coupling limit, where A 
is large. Thermodynamic quantities like the free energy in particular 
were computed in both theories and they do not agree with each other 
for this reason. 

7 Discussion 

The phase transition in the matrix model is a peculiar one. Closer 
analysis reveals that it is of first order. It is easy to see from (27) 
that the expectation value of the Polyakov loop operator is zero in the 
symmetric phase but is non-vanishing in the high temperature phase 
and approaches a non-zero value there even as the critical point is 
approached from above. However, at this critical point, there is no co- 
existence region where one phase is meta-stable and the other is stable, 
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as is normally the case for first order phase transitions. To see this, 
note that the potential the phase transition occurs where the potential 
becomes unstable to perturbations. Normally, in a first order behavior, 
there are two competing vacua, both of which are perturbatively stable 
in the transition region and as parameters are varied one or the other 
gets lower free energy and is preferred. Then they cross there is a phase 
transition. In the present case, the perturbative instability occurs at 
the same place as the phase transition. 

There is a question as to whether this is an artifact of the ap- 
proximation that we have done here, i.e. if we expanded the effective 
action for eigenvalues to higher loop order the phase transition might 
be more conventional. 

Indeed, in other unitary matrix models applied to non-Abelian 
Coulomb gases [34, 35], where the eigenvalues live on a higher dimen- 
sional space (a D = 1 or even D > 1 unitary matrix model) there is a 
coexistence region. 

A similar problem afflicts weakly coupled four-dimensional super- 
symmetric Yang-Mills theory [15]. In that case, Witten argued using 
AdS/CFT that, in the strong coupling limit of planar Yang-Mills the- 
ory, there should be a de-confining phase transition, identified with 
the Hawking-Page phase transition of supergravity on asymptotically 
anti-de Sitter space [36]. Hawking-Page is a normal first order phase 
transition where there is the possibility of metastable phases. At weak 
coupling, the analysis looks very similar to what we have done here 
for the matrix model and the phase transition has the same nature. 
Aharony et.al. have conjectured that the effect of higher loop cor- 
rections (3-loops) to the effective action in Yang-Mills theory would 
indeed change the phase transition to a more conventional one. 

Finally, in the matrix model that we have analyzed, there is the 
question of whether the phase transition that we have found has any- 
thing to do with collapse to black holes. The subject of black holes on 
pp-wave backgrounds is a murky one which is presently being sorted 
out. In any case, the physics described by classical gravity should 
appear at a strong coupling limit, rather than the weak coupling limit 
that we have analyzed. It is tempting to conjecture that, if there were 
black holes, collapse to black holes at finite temperature is what our 
phase transition would describe if it persists at strong coupling. 

There is the further question of whether our phase transition is 
related to the Hagedorn behavior that is seen in the string spectrum 
on pp-wave backgrounds [37, 38, 39, 40, 41, 42, 43]. This behavior 
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of course occurs in weakly string theory which is a limit of the M- 
theory. It there were such a relationship, it would be interesting to ask 
whether it is related to the formation of black holes in ten dimensional 
supergravity near a pp-wave background [44, 45]. 

Of course of the conjectured relationship with the 5-brane is valid, 
then, indirectly, the Hagedorn phase transition of little string theory is 
related to horizon formation in the 5-brane geometry. In that case, our 
phase transition and its thermodynamics indeed describes the black 
hole, again in a limit which is far away from previous analysis of such 
objects. Note that we do not analyze the relative stability of mem- 
branes and five-branes. This question has been addressed in recent 
interesting papers[46, 47] 

It is also interesting to ask whether the phase transition that we 
have identified is related to recent work which studied phase transi- 
tions for statistical systems of random walks on discrete groups like 
the permutation group. [48, 49] 
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